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ABSTRACT 

Abundances of C, N, and O are determined in four bright red giants that span the known abundance range 
for light (Na and Al) and i-process (Zr and La) elements in the globular cluster NGC 185L The abundance 
sum C+N+O exhibits a range of 0.6 dex, a factor of 4, in contrast to other clusters in wh ich no significant 
C+N+O spread is found. Such an abundance range offers support for the lCassisi et al.l ( l2008h scenario in which 
the double subgiant branch populations are coeval but with different mixtures of C+N+O abundances. Further, 
the Na, Al, Zr, and La abundances are correlated with C+N+O, and therefore, NGC 1851 is the first cluster to 
provide strong support for the scenario in which AGB stars are responsible for the globular cluster light element 
abundance variations. 

Subject headings: Galaxy: Globular Clusters: Individual: NGC 1851, Galaxy: Abundances, Stars: Abundances 
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1. INTRODUCTION 

Globular clusters have long been regarded as simple stel- 
lar populations which may be described by a single age, he- 
lium abundance (Y), metallicity (Z), and initial mass function 
(jkenzini & Buzzoni 1986). For many years these simple stel- 
lar populations have played a prominent role in astronomy 
by providing a lower limit to the age of the Universe and for 
testing the predictions of stell ar evolution and ste llar nucle- 
osynthesis (e.g., see review by lGratton et al.ll2004l) . 

However, a revolution in the field of globular cluster re- 
search is underway, prompted by the recent discoveries of 
complex structure on the main sequence, subgiant branch 
(SGB), red giant branch (RGB), and/or horizontal branch 
(HB) within some Galactic and extra-Galactic globular clus- 
ters indicating that these simple stellar populations in fact 
contain discrete, mult i ple populations (e. g., Bedin et al. 200' 
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D'Antona et alj|2005l: iPiotto et al.ll2005LSonima et al...2005l 
Mackev et alJ l2008l) . At present, these multiple populations 



can be best explained in terms of distinct ages and/or com- 
positions, although the sequence of eve nts leading to their 
formation remains largely unexplained dRenzinil i2(j08l) . In 
some cases, extremely large heUum abundances, up to F = 
0.40, are required to explain the vario us populations (e.g., 
ID' Antona & Cal o? '2004?: 'Norris' '200?). While theoretical 
mode ls have strug gled to account foi^ such high abundances 
(e.g.. lKarakas"etalJi2006ai: iBekki & Chiball2007|). some suc- 
cessful efforts are begin ning to appear (e.g., iD'Ercole et alJ 
120081; iPumo et aLll2008l) . The intrigue surrounding the mul- 
tiple populations in globular clusters has been enhanced by 
the speculation that clusters with extended HBs, generally 
the most massive clusters, are the remnants of dw arf galax- 
ies from which the halo of the Galaxy was built dLee et alj 

' Based on observations made with the Magellan Clay Telescope at Las 
Campanas Observatory. 
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l2007h . 

The globular cluste r NGC 18 5 1 has a bimodal HB and dis- 
plays multiple SGBs (iMilone et al. 2008, hereafter M08), but 
no evidence for multiple main sequences. Based on the width 
of the main sequence and RGB, M08 suggested that the max- 
imum possible abundance variation is A[Fe/H] = 0.1 dex or 
Ay = 0.026. From an analysis of eight bright RGB stars using 
high resolution spectroscopy. lYong & Grunda hl (2008, here- 
after YG08) showed that the dispersion in [Fe/H] was less 
than 0.1 dex, and that NGC 1851 displays the usual star-to- 
star abundance variations for the light elements O, Na, and 
Al; the exact mechanism for these abundance patterns, seen 
in all globular clusters, remains unknown (e.g.. lGratton et al] 
lIOOl . However, YG08 also showed that NGC 1851 harbors a 
star-to-star abundance variation for the .s-process elements Zr 
and La and that the abundances of these elements were corre- 
lated with the light element abundances. Furthermore, within 
the small sample there was a hint that the abundances of the s- 
process elements was bimodal, which suggests that the RGB 
may consist of two populations with distinct chemical com- 
positions. Therefore, as well as showing unusual photometric 
properties, NGC 1851 also displays peculiar chemical abun- 
dance patterns not seen in other clusters. 

Several explanations for the double SGB in NGC 1851 have 
been offered. If age is the sole parameter, M08 suggested 
a 1 Gyr difference with the fainter subgiant branch (fSGB) 
population being older M08 also considered a combination 
of increasing [Fe/H] by 0.2 dex and helium from Y = 0.247 to 
0.30, though this possibility was excluded by the magnitudes 
of the blue HB s t ars. A n alternative explanation was proposed 
by lCassisi et al.l (12008 ) in which the two SGB populations are 
coeval, but with the fSGB population having a total C+N+O 
abundance increased by a factor of 2. In this scenario, the 
two populations may have comparable He abundance (at most 
Ay = 0.032), in contrast to the large He variations inferred in 
other clusters which display multiple populations. We have 
obtained spectra with the goal of m easuring the abundance 
sum C+N+O to test the lCassisi et alJ (120081) scenario. 

2. OBSERVATIONS AND ABUNDANCE ANALYSIS 

We selected four bright RGB stars from YG08, two Na-, 
A1-, Zr-, and La-normal targets as well as two Na-, A1-, Zr- 
, and La-rich objects. In Figure [T] we show the locations of 
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Fig. 1. — The color-magnitude diagram for NGC 1851 using the 
IGru ndahl et al.' (1999) photometry. The large plus symbols represent the 
bright giants in this study. The filled squares represent the additional four 
stars in YG08. 

our targets in a color-magnitude diagram. If the double SGB 
populations are due to distinct chemical compositions, and 
if these chemical compositions are also present on the RGB, 
then our small, but carefully selected, sample should allow us 
to identify any C-i-N-i-O abundance variation, if present. 

The stars were observed with the high resolution spectro- 
graph MIKE (Bernstein et al. 2003) on the Magellan Clay 
Telescope on 2007 December 23. We used the 0.7" slit which 
provided a spectral resolution of R ^ 37,000 in the blue CCD 
(3350-5000A) and R ~ 30,000 in the red CCD (4900-9500A) 
as measured from the ThAr lines. Our total exposure times 
per star ranged from 40 to 80 minutes resulting in signal-to- 
noise ratios (S/N) per four pixel resolution element of 70-100 
at4320A and 210-260 at 6600A. The data were redu ced using 
the same procedure described bv lYong et al.l (l2006h . 

Although the new data are of lower spectral resolution {R ~ 
30,000 vs. R ~ 55,000), the S/N and wavelength coverage are 
superior to the spectra analyzed by YG08. Using these new 
spectra, we rederived the stellar parameters using the same 
tools and techniques described in YG08 and found that the 
revised stellar parameters are in good agreement with the stel- 
lar parameters previously obtained. The new dispersion in Fe 
abundances for these four stars is cr[Fe/H] = 0.03 dex com- 
pared with 0.11 dex for the YG08 analysis. This decrease is 
due to the increased wavelength coverage and S/N which al- 
lows us to more accurately measure the equivalent widths of 
a larger number of Fe lines. 

The abundances of C and N were derived by comparing 
observed spectra with synthetic spectra generated using the 
LTE line analy sis and spectrum synthesis program MOOG 
( ISnedenlll973b . For C, we analyzed lines from the (0,0) 
and (1,1) bands of the A-X electronic transition of the 
CH m olecule near 431 OA using the line list from iPlez et alj 
(120081) . For N, we analyzed lines from the (2,0) band of the 
A-X electronic transition of the C N molecule near 8000A 
usi ng the line list from Reddv et al. ( 2002) which is based on 
the Ide Lavernv & Gustafsson (1998 ) list. In both cases, we 
adjusted the abundances of C or N until the synthetic spectra 



matched the observed spectra (see Figures |2] and |3]l. We red- 
erived the O abundances again from the 6300A [OI] line using 
the new spectra, the revised C and N abundances, and the re- 
vised stellar parameters. Since the abundances of C, N, and 
O are coupled, we iterated until self-consistent abundances 
were obtained. We also updated the abundances for Na, Al, 
Zr and La using the new spectra and revised stellar parameters 
(a complete abundance analysis of our data will be presented 
elsewhere). We present the stellar parameters and abundances 
in Tabled 
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Fig. 2. — Observed spectra (circles) and synthetic spectra (lines) for C 
(upper), N (middle), and O (lower) for Star 395. The synthetic spectra show 
the best fit (thick black line) and unsatisfactory fits (thin red and green lines) 
A(C) ± 0.30 dex, A(N) ± 0.20 dex, and A(0) ± 0.20 dex. 



3. RESULTS 

As mentioned, our sample was selected to span the full 
range of the known star-to-star abundance variation of Na, Al, 
Zr, and La. Our new analysis confirms the large abundance 
spread for these elements as well as the correlations between 
the light (Na and Al) and heavy (Zr and La) elements estab- 
lished by YG08. The range in Na, Al, Zr, and La abundances 
far exceeds the measurement uncertainties of 0.06, 0.05, 0. 12, 
and 0.08 dex respectively, though we note that the abundance 
amplitudes for these elements are slightly smaller than those 
reported by YG08. 

Adopting a solar a bundance log e(Fe) = 7.50 
(iGrevesse & Sauvail Il998l) . we find a mean cluster abun- 
dance of [Fe/H] = -1.22 (tr = 0.03). The revised cluster Fe 
abundance is in good agreement with the YG08 value and 
previous estimates in the literature, but the abundance dis- 
persion is considerably smaller for the reasons noted above. 
The abundance range for Fe may be entirely attributable to 
observational uncertainties (0.07 dex). 
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TABLE 1 

Stellar parameters and abundances 



Star" 


V 


Teff (K) 


\ogg (cgs) 


IV (km s 


') A(Fe) 


A(C) 


A(N) 


A(0) 


A(CNO) 


A(Na) 


A(A1) 


A(Zr) 


A(La) 


003 


13.68 


4125 


0.7 


1.86 


6.28 


6.46 


7.60 


7.75 


8.00 


5.25 


5.63 


1.51 


0.16 


095 


13.58 


4125 


0.7 


1.92 


6.32 


6.51 


7.60 


7.60 


7.92 


5.25 


5.69 


1.64 


0.22 


329 


13.46 


3925 


0.0 


1.92 


6.25 


6.81 


8.00 


7.88 


8.26 


5.50 


5.79 


1.87 


0.42 


395 


13.48 


4025 


0.5 


2.01 


6.29 


6.41 


8.45 


7.38 


8.49 


5.74 


6.04 


1.85 


0.53 



' Star names taken from lStetsoj fT98lh . 
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Fig. 3. — Same as Figure |2]but for Star 95. 



For C and N, we find a large range in abundances. The 
0.40 dex range in log e(C) and the 0.85 dex range in log e(N) 
greatly exceed the measurement uncertainties (0.04 dex and 
0.11 dex respectively). The abundance amplitude for N, 
A log e(N) = 0.85, is smaller than that seen in other clus- 
ters, e.g., RGB bump stars in NGC 6752 and main sequence 
stars in Ml 3 show 2.0 dex amplitudes (iBriley et alJ 120021 : 
lYongetalJl2008al) . This may be due to the small sample 
and/or advanced evolutionary status where extensive CN- 
cycling may reduce the amplitude of the N abundance disper- 
sion. Curiously, we do not detect any obvious anticorrelation 
between C and N as seen in every well studied Galactic glob- 
ular cluster We identify correlations between the abundances 
of N and Na as well as N and Al, but the corresponding anti- 
correlations between C and Na and C and Al are not obvious. 
Nevertheless, we find correlations between the abundances of 
N and Zr as well as N and La as expected given the known 
behavior of N and Al and the correlations between Al and Zr 
and Al and La found YG08. 

Our first significant result is that the abundance sum 
Ch-Nh-0 shows a large variation. Within our small but bi- 
ased sample, the amplitude of the abundance variation is 
ZU(Ch-Nh-O) = 0.57 dex. The 0.57 dex spread far exceeds 
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Fig. 4. — The abundances of Na, Al, Zr, and La vs. the abundance sum 
C+N+O. A representative error bar and a linear least squ ares fit to the data 
(including form al slope and erro r) are shown. Data for M4 jSmith et al.l200^ 
and NGC 6712 jYong et aU2008b ^ ai'e included for compaiison. 

the estimated uncertainty of 0.14 dex. Bright giants in M4 
and NGC 6712, globular clusters of comparable metallicity 
to NGC 1851, show Ch-Nh-O abundance ampHtudes o f 0.22 
dex an d 0.35 dex respectively (Smith et al. 2005; Yon g et al] 
I2008bl) . Unevolved stai's in NGC 6397, NGC 6752 and 47 
Tuc show Ch-Nh-O abundances that ma y be regarded as con - 
stant within the ^0.3 dex uncertainties (ICarretta et al.ll2005h . 

Our second main result is that the abundances of Na, Al, 
Zr, and La show a positive correlation with the abundance 
sum Ch-Nh-O (see FigurelDi. In this figure, we overplot results 
for M4 and NGC 6712 and find that for the available data, 
there are no positive correlations between A(X) and C+N+O 
in these other clusters. Finally, our small sample does not 
exclude the possibility that the C+N+O abundance may be 
bimodal within this cluster 

4. DISCUSSION 

We seek to test the lCassisi et alJ (|2008|) prediction in which 
the double SGB in NGC 1851 is due to two coeval popu- 
lations whose chemical compositions for C+N+O differ by 
a factor of two. Note that we have not directly observed 
stars on the SGB. Instead, we have observed a sample of 
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bright RGB stars in which we have found a large spread in 
the C+N+0 abundances. If the SGBs comprise populations 
with different C+N+0 abundances, then we assume that these 
populations with distinct compositions are also present on the 
main sequence and R GB. Stromgren photometry obtained by 
iGrundahl et alJ (11999) reveals that the bimodality seen on the 
SGB can also be found on the RGB. Therefore, the large 
spread in C+N+O abunda nces in RGB stars o ffers support, 
albeit not definitive, for the lCassisi et al.l (120081) scenario. We 
tentatively confirm that the double SGB in NGC 1851 is due 
to two coeval populations with C+N+O abundances differing 
by a factor of two. 

Based on the relative numbers of stars, YG08 speculated 
that the brighter subgiant branch (bSGB) stars were related to 
the RGB stars with "normal" Zr and La abundances and that 
these stars populate the red HB. YG08 also suggested that the 
fSGB stars are related to the RGB stars with Zr and La en- 
hancements and that these stars populate the blue HB. If the 
fSGB is older than the bSGB as M08 suggest, then the YG08 
speculation would require the Zr-rich, La-rich stars to be older 
than the Zr-normal, La-normal stars. In this scenario, the se- 
quence of events leading to cluster formation would be com- 
plicated and require the AGB (or other) stars to preferentially 
pollute the older Zr-rich La-rich population but without pol- 
luting the gas from which the younger Zr-normal La-normal 
population formed. However, if the C+N+O abundance of the 
i-process rich stars is larger than that of the i-process nor- 
mal ones, the fSGB may be coeval to the bSGB, or even a 
bit younger than the bSGB, as shown in Cassisi et al. (2008). 
Due to the observational errors, and to the fact that the age 
difference we are looking for does not exceed some lOOMyr, 
it may be very difficult or impossible to constrain the precise 
mass range of the progenitors of the i-process rich stars, but 
a coherent picture of the chemical evolution history is emerg- 
ing. 

The ICassisi etal.1 (l2008l) scenario has important impli- 
cations for other cl usters which display multiple SGBs. 
IVillanova et al.l (l2007i) studied the SGB region of uo Centauri, 
the most massive globular cluster which is suspected of har- 
boring up to five populations with discrete ages and/or com- 
positions. From the large range in the magnitudes of the SGB 
populations, iVillanova et al.l (120071) inferred a large range in 
ages. They found that there was an extremely old, but metal- 
rich population which greatly complicates any formation sce- 
nario. Based on our NGC 1851 results, the inferred ages of 
SGB stars in u) Cen requires detailed knowledge of the rela- 
tive abundances of C+N+O. We speculate that the old metal- 
rich population (Villanova et al.'s "SGB group D") has a large 
excess of [(C+N+0)/Fe] relative to the more metal-poor pop- 
ulations and therefore a younger relative age than currently 
inferred. Such compositions would be consistent with our un- 
derstanding of the composition of the ejecta from AGB stars 
dKarakas & Lattanzio 2007 ) and the pro minent role of AGB 
stars in the chemical evolution of uo Cen (iNorris & Da Costal 
[1991 . 

AGB stars have long been suspected as being responsible, 
in part or in whole, for the star-to-star abundance variation for 
the light elements (C, N, O, F, Na, Mg, and Al) found in every 
well studied Galactic globular cluster, including NGC 1851. 
However, various theoretical yields and chemical evolution 
models involving AGB stars ha ve yet to provide a satisfac- 
tory solution (e.g. jFenner et al .''2004). AGB stars may be ex- 
pected to produce a substantial increase in the C+N+O abun- 
dance as they enhance Na and Al and deplete O and Mg. AGB 



stars should also produce i-process elements, albeit via differ- 
ent processes than for the O, Na, Mg, and Al nucleosynthesis. 
If AGB stars are responsible for the light element abundance 
variations in globular clusters, we may expect a large range 
in C+N+O abundances along with correlations between (1) 
C+N+O and Na, (2) C+N+O and Al, and (3) C+N+O and s- 
process elements. On the other hand, these proposed abun- 
dance trends are dependent on the number of third dredge 
up episodes and therefore on the adopted convection effi- 
ciency (e.g., Ventura & D'Antona 2005a; Renzini 2008) and 
mass loss (e.g., V entur a & D'Anton a 2005b). Further uncer- 
tainty in the yields c omes from the uncertainties in the nuclear 
reaction rat es (e.g., Lusaro et al."2004'; 'Vent ura & D'Antonal 



2006; Karakas et al. 



2006b; Izzard et al. 20071). Although the 



brightest AGB stars in the LMC and SMC can become C-rich, 
it is possible that nucleosynthesis in the most massive AGBs 
is scarcely affected by the third dredge up, with only a minor 
increase in the total C+N+O abundanc e of the second gener- 
ation stars (iD'Antona & Ventural2008l) . Nevertheless, the ab- 
sence of the signatures listed above has forced astronomers 
to consider alternative sources such as massive stars (e.g., 
iDecressiiTet al. 2007). 

The correlations between C+N+O and Na, Al, Zr, and La 
presented here offer strong support for the scenario in which 
AGB stars are responsible for the light element abundance 
variations. Given the peculiar photometric and chemical prop- 
erties of NGC 1851, it is surprising that this cluster offers the 
strongest support for the globular cluster AGB pollution sce- 
nario. It would be ironic if NGC 1851 is the only globular 
cluster for which the AGB pollution scenario is applicable. 
An alternative interpretation is that the AGB pollution sce- 
nario is applicable to all globular clusters. For NGC 1851, 
the masses of the polluting AGB stars was lower such that the 
products of third dredge-up, increased C+N+O and s-process 
elements, are clearly evident. 

We note that abund ance spreads in A('CNO) are likely 
present in uo Centauri dNorris & Da Costal 1 1995h . although 
they were concerned with possible problems with the N mea- 
surements. Indeed, NGC 1851 bears a striking resemblance 
to u) Centauri. Both clusters possess multiple SGBs, large 
abundance spreads for ^-process elements and likely spreads 
in A(CNO), and a scatter in the Strom gren m\ index , tradi- 
tionally used as a metallicity indicator ( Grundahl et al.iri999,) . 
We reiterate the comments by YG08 that NGC 1851 is a 
"bridge" between uj Centauri (variation of all elements) and 
NGC 6752-like clusters (constant Fe and C+N+O but large 
variations of light element s C-Al). 

Finally, recent work by ISalaris et al.l (l2008l) also suggests 
a relationship between the SGB and HB populations in NGC 
1851. Specifically, the bSGB stars evolve onto the red HB 
while the fSGB stars populate the red and blue HB. If this 
is correct, we would therefore expect stars on the red HB to 
be primarily Zr-normal and La-normal. NGC 1851 may pro- 
vide critical information on HB morphology and the second 
parameter effect and additional studies of this cluster are of 
great interest. 
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